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ABSTRACT 

The investigation was an approximate method for the analysis of rigid framed medium rise 
building subjected to aerodynamic loadings. The main aim of the research work is to develop 
appropriate equations and formulae that will enable the use of the building variable dimension 
(i.e. such as Li, L 2 , h etc.) and either table calculators or other related computer spread sheet to 
determine the axial, shearing forces and their corresponding bending moments on the columns 
and beams; at any level of a rigid unbraced frame subjected mainly to aerodynamic loadings. It 
assumed that at the middle of columns and beams, are located points of contra flexures. The 
axial forces caused by the effects of the lateral loadings in the columns are proporţional to their 
distances from the center of gravity of the frame. For low and medium rise buildings, the 
columns are assumed to be of the same cross-sectional area. For buildings lower than five 
storeys, wind loadings are negligible since the floor slab or diaphragm is assumed to be 
infinitely rigid but for those higher than five storeys, the effects of wind loadings are 
considered important for the analysis of the structural system. At heights greater than lOm, the 
wind speed is assumed to increase approximately, linearly with height. Therefore, a coefficient 
k, was incorporated to take into consideration the increase of the basic wind speed with height, 
local site conditions, building category, directionality and height coefficient respecţi vely. 
Useful formulae and results obtained show the approaches are appropriate methods that can be 
used to obtain approximate preliminary sections suitable for analysis and design of beams and 
columns for low and medium rise multi-storey building subjected to wind gust anywhere in the 
World. 
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INTRODUCTION 

Rigid Frame Multi-Storey Buildings and Methods of Analysis 

Rigid-frame high-rise buildings typically consist of a conglomeration of parallel or orthogonally arranged structural 
elements such as columns and beams or girder with moment resisting joints. Resistance to horizontal loadings is 
usually provided by the flexural resistance of the slab diaphragms, columns, girders and the monolithic joints. The 
continuity of the frame also contributes to the resistance from the vertical loadings such as dead and superimposed 
live loads by reducing the bending moment in the beams or girders. Simplified practicai examples of this method of 
design for columns have been elaborately discussed by many authors (Oyenuga 2011, Mosley et a/.;2007) 


According to Smith and Coull (1991), rigid frame is economical only for buildings up to about 25 stories; above 25 
stories, the relatively high lateral flexibility of the frame will require very uneconomically large members for the 
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drift at higher floors to be effectively controlled. Gravity loading is also resisted by the rigid-frame action (Oyenuga 
2011). Negative moments are induced in the girders adjacent to the column causing the midspans positive moments 
to be significantly reduced when compared with midspans moments from equivalent spân of a simply supported 
beam. Regardless of this advantage, in structures where gravity loads dictate the design concept, economies that 
arise from this effect tend to be offset by the high cost required to actualize effectively monolithic rigid joints. 

According to Grider et al; (1999), a number of methods have been developed over the years for the analysis of 
continuous beams and frames. The so-called classical methods — such as application of the theorem of three 
moments, the method of least work, and the general method of consistent deformation have proved useful mainly in 
the analysis of continuous beams having few spans or of very simple frames. For the more complicated cases usually 
met in practice, such methods prove to be exceedingly tedious, and alternative approaches are preferred. For many 
years the closely related methods of slope deflection and moment distribution provided the basic analytical tools for 
the analysis of indeterminate beams and frames. In offices with access to high-speed digital computers, these have 
been supplanted largely by matrix methods of analysis. Where computer facilities are not available, moment 
distribution is still the most common method. Approximate methods of analysis, based either on an assumed shape 
of the deformed structure or on moment coefficients, pro vide a means for rapid estimation of internai forces and 
moments. Such estimates are useful in preliminary design and in checking more exact Solutions and in structures of 
minor importance may serve as the basis for final design. 

The method of slope deflection entails writing two equations for each member of a continuous frame, one at each 
end, expressing the end moment as the sum of four contributions: (1) the restraining moment associated with an 
assumed fixed-end condition for the loaded spân, (2) the moment associated with rotation of the tangent to the 
elastic curve at the near end of the member; the moment associated with rotation of the tangent at the far end of the 
member, and the moment associated with translation of one end of the member with respect to the other. These 
equations are related through application of requirements of equilibrium and compatibility at the joints. A set of 
simultaneous, linear algebraic equations results for the entire structure, in which the structural displacements are 
unknowns. Solution for these displacements permits the calculation of all internai forces and moments. This method 
is well suited to solving continuous beams, provided there are not very many spans. Its usefulness is extended 
through modifications that take advantage of symmetry and anti-symmetry and of hinge-end support conditions 
where they exist. However, for multistory and multi-bay frames in which there are a large number of members and 
joints, and which will, in general, involve translation as well as rotation of these joints, the effort required to solve 
the correspondingly large number of simultaneous equations is prohibitive (Grider et al.; 1999). 


Use of matrix theory makes it possible to reduce the detailed numerical operations required in the analysis of an 
indeterminate structure to systematic proces ses of matrix manipulation which can be performed automatically and 
rapidly by computer. Such methods permit the rapid solution of problems involving large numbers of unknowns. As 
a consequence, less reliance is placed on special techniques limited to certain types of problems; powerful methods 
of general applicability have emerged, such as the matrix displacement method. Account can be taken of such 
factor s as rotational restraint provided by members perpendicular to the plane of a frame. A large number of 
alternative loadings may be considered. Provided that computer facilities are available, highly precise analyses are 
possible at lower cost than for approximate analyses previously employed. 


In spiţe of the development of refined methods for the analysis of beams and frames, increasing attention is being 
paid to various approximate methods of analysis. There are several reasons for this. Prior to performing a complete 
analysis of an indeterminate structure, it is necessary to estimate the proportions of its members in order to know 
their relative stiffness upon which the analysis depends. These dimensions can be obtained using approximate 
analysis. Also, even with the availability of computers, most engineers find it desirable to make a rough check of 
results using approximate means to detect gross errors. Further, for structures of minor importance, it is often 
satisfactory to design on the basis of results obtained by rough calculation. Provided that points of inflection 
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(locations in members at which the bending moment is zero and there is a revers al of curvature of the elastic curve) 
can be located accurately, the stres s resultants for a framed structure can usually be found on the basis of static 
equilibrium alone. Each portion of the structure must be in equilibrium under the application of its externai loads 
and the internai stres s resultants. 

Mo st of the methods mentioned above involve many hours of rigorous analysis and if care is not taken, inherent 
mistakes and errors could creep into the computaţional process. Even the use of computer methods and soft ware’s 
on most occasions require the entry of large volumes of input and output data or the development of complex 
algorithms for Solutions of multi-storey framed buildings. But for reliable Solutions to be obtained, preliminary 
sizing of main structural members must be determined by approximate methods and entered as input. Therefore, the 
main aim of this research work is to develop appropriate equations and formulae that will enable the use of the 
building variable dimension (i.e. such as L ls L 2 , h etc.) and either table calculators or other related computer spread 
sheet to determine the axial, shearing forces and their corresponding bending moments on the columns and beams; 
at any level of a rigid unbraced frame subjected mainly to aerodynamic loadings. These formulae are mainly 
applicable to the full or preliminary assessment of structural members for: 1) For full assessment of medium rise 
multi-storey buildings that satisfies height / width ratio less than six, overall height less than 50m and subjected 
mainly to aerodynamic loadings. 2.) They are also useful for the preliminary assessment of buildings that have 
dynamic response with minimum frequency (f < 1.0 Hz) or fundamental natural period exceeding one second (T > 1 
sec) (Onundi et al.; 2012). 

Wind and Wind Loadings 

Taranath (2005), wind is a term used to describe horizontal motion of air. Motion in a vertical direction is called a 
current. Winds are produced by differences in atmospheric pressure that are primarily attributable to differences in 
temperature. These differences are caused largely by unequal distribution of heat from the sun, and the difference in 
thermal properties of land and ocean surfaces. When temperatures of adjacent regions become unequal, the warmer, 
lighter air rises and flows over the colder, heavier air. Winds initiated in this way are modified by rotation of the 
earth. Near the equator, the lower atmosphere is warmed by the sun’s heat. Forces due to wind are generated on the 
exterior of the building based on its height, local ground surface roughness (hills, trees, and other buildings) and the 
square of the wind velocity. The weight of the building, unlike in earthquake design, has little effect on wind forces, 
but is helpful in resisting uplift forces. Unless the structure has large openings, all the wind forces are applied to the 
exterior surfaces of the building. This is in contrast to earthquake forces where both exterior and interior walls are 
loaded proportionally to their weight. Wind pressures act inward on the windward side of a building and outward on 
most other sides and most roof surfaces. Special concentrations of outward force, due to aerodynamic lift, occur at 
building corners and roof edges, particularly so at overhangs or eaves of the building. The overall structure is 
designed for the sum of all lateral and uplift pressures and the individual parts to resist the outward and inward 
pressure concentrations. They must be connected to supporting members to form a continuous load path. 

According to Kijewski and Kareem (2001), Mendis et tf/.;(2007) and Joseph and William (2006), many aspects 
involved in the estimation of wind loads are held in common by the various internaţional codes and standards. 
Instead of commenting on them repeatedly, they have been highlighted in their work briefly that:: 

a) All the standards, subdivided the global terrain into 3 to 5 categories depending on how they affect the 
wind characteristics at that location. 

b) The design wind speed, associated with one or a range of mean recurrence intervals, used in analysis 
by each of the codes is typically the product of the basic wind speed and factors to account for the 
geographic location, topographical effects, building size and surface roughness, etc. 

c) Wind gustiness introduces dynamic load effects which the codes and standards account for by 
factoring up the mean loads by a gust factor. Both time and spaţial averaging play an important role in 
the development of gust factors. For a very small size structure, a short duration gust, which 
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completely engulfs the structure, e.g. a 3 second gust, may be adequate to account for the effects of 
gustiness, in which case the gust factor is unity. 

d) On the other hand, if the wind-averaging interval is higher, e.g. 10 minutes or more, the averaged wind 
exhibits less fluctuation, and accordingly the gust factor is greater than unity. This departure from 
unity is affected not only by the averaging interval, but also by the site terrain and the size and 
dynamic characteristics of the structure. 

e) Furthermore, it should also be noted that while all of the standards reference their wind speed at 10 m 
above ground in a flat, open exposure, each uses gusts of different duration. The British and Canadian 
standards use the mean hourly wind speed in design, while the European Pre standard, the China 
National Standard and the AIJ Recommendations all use a 10-minute mean wind velocity. The 
ASCE7-95 standard references a 3 second gust, as does the Australian Standard, though, in the latter 
case. This wind is later converted to a mean hourly wind for subsequent calculations of dynamic 
pressure and the gust factor. As a result, for any adequate comparison amongst standards, there must 
be proper adjustments of the reference velocity. 

Determination of the aerodynamic loading on the Building 

The dynamic wind load, F along the height, L of a high rise building is given as 

Fi = Pe A (1) 


Where; 

Pe = the pressure acting on the externai surface of a building 

P e = q s C pe C a C r (BS 6399 (2004) (la) 

C pe = the externai pressure coefficient for the building surface given in clauses 2.4 and 2.5 
of the Code; C a = the size effect factor for externai pressures defined in clause 2. 1.3.4. 
q s = the dynamic pressure from clause 2.1.2 and q s given as 0.6 13 V e 2 

V e = the effective average basic wind speed in m / sec (Soboyejo, 1971, Onundi et al.; 2009 and Onundi, 2010) 
and 

A = the loaded area 

••• Fi = q s AC pe C a C r (lb) 

The dynamic augmentation factor, (C r ) recommended by the Clause 1.6.1, Figure 3; BS 6399 (2004) for the high 
rise building. The lateral loadings acting on a structure may be caused by either wind pressures, earthquakes (i.e. 
seismic forces) or earth retained by the walls of the structure. The method described in this study assumes that the 
vertical loadings are computed separately and will be combined with the results obtained from the effects of the 
lateral loadings. 


Analysis of Medium Rise Buildings 

For the design of low to medium high rise buildings, the elastic method of analysis and the folio wing approximation 
methods of analyses may be assumed (Mosley et al; 2011). 

i. At the middle of columns and beams, are located points of contra flexures. 

ii. The axial forces caused by the effects of the lateral loadings in the columns are proporţional to their 
distances from the center of gravity of the frame. For low and medium rise buildings, the columns are 
assumed to be of the same cross-sectional area. 

iii. For buildings lower than five storeys, wind loadings are negligible since the floor slab or diaphragm is 
assumed to be infinitely rigid but for those higher than five storeys, the effects of wind loadings are 
considered important for the analysis of the structural system. At heights greater than lOm, the wind speed 
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is assumed to increase approximately, linearly with height. Therefore, a coefficient k, is incorporated to 
take into consideration the increase of the basic wind speed with height, local site conditions, building 
category, directionality and height coefficient respectively as shown in Figure 1.0(a). 


By using theprinciple of similar triangle, it becomes possible for intermediate values, k x between k max and k^n to be 
determined from the following eq 




K x _ K ma x 
L-x L k 


Therefore, 

K x =K max ^ (2) 

If L x = 0; 
kx — k min — k —0 

and W min = k + W = W. (2a) 


This will correspond to the basic wind speed W, at lOm height 
for the site. 

If 0 < L x < L k , 
k x = k int 

and Wint — k x + W m m — ki n ţ + ~WrniTr (2b) 

This will correspond to the wind speed intermediate 
between lOm and the maximum height L k , of the building. 

If L x = L k , 
k x = k max and 

Wmax — k x + W m j n — k max + W)^. (2c) 

This will correspond to the maximum wind speed 
corresponding to the maximum height of the building. 


Figure 1.0 (a) : Linear Variation of Wind Speed along the Building 

a.) At the Top Floors 





To determine the nodal force Fi at any storey height i, the following 
Procedures are suggested. 

Let the force at top floor be represented by F n and W max be the 
expected maximum wind pressure at the building top floor Figure 1.0(b) 
The corresponding wind pressure at the building intermediate 
storey heights between the top floor and lOm height measured 
from ground level and is represented by W int and the nodal forces 
are also determined by F n _ m . 

Figure 1.0(b) : Linear Variation of Wind Speed above lOm Height 
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Let W max = (K 

max + W min )b 
Win, = (K int + W mi „)b 

But 

K X| _ K max _ jh 

x ’~^ max 

K xj = K j =iK max =^K max (2d) 

n = total number of storeys in the mulţi- storey building 

j = total number of storeys above the datum level (i.e. in this case, the datum is lOm) 
b = Width or frame spacing 

Therefore, 

Wj = (w min + iK max )b 
or 

Wj = (W min + 2±iK max )b (2e) 

Wj = Wind speed at j th storey level 

Alternaţi vely, for all the floors, when n > j > lOm 
Kn = = K min + K n = K min + ^K max and 

W„ = W max = (W min + ^ K max )b = 0.0.1042 b W bas gj) (2f) 

W bas = Basic wind speed at the building site. 


For the floors lower than lOm, when j < lOm 

K' = K n 


+ K n = 0 + K n 


Wj = W n = bW min 

Therefore, 

Figure 2.0(a, b,c and d), shows that at the top floor, the nodal force is 


.(2g) 


F„=w!L( 


Kmax + Kj n t+ 2K min 


) — W^-(K n + K m + 2K n 


J = K'wŞ = W n ^ 


Whereas, the corresponding wind pressure at the building intermediate storey heights (i.e. floors) between the top 
floor and the five storey height have the nodal forces represented by F n+ i_ m 

h 2 /K n _ m + K m + 2K lr , in \ h 2 . h 2 h 2 h 2 


h /r 

F n m = W — - 
n “ m 2 V 

yj ^ K m + 2K m j n ^ 


w 


h 2 


) = W-(K n _ m + K m + 2K r 

) 


J = k' — = K' — = Wi — 

ny lv n-m vv A lv n-m vv min A vv j A 


fK„ ^ 

v lv n— m 1 lv m 1 


KmW 


h 2 


K' 

lv m v v min 


h 2 


W r 


h 2 


4 

....( 3 ) 


Where 

K n = the maximum variable wind pressure coefficient corresponding to the building top floor 
(e.g. if n=10, hence K n = Ki 0 ) 

K n . m = the intermediate variable wind pressure coefficient corresponding to the building 

intermediate floors counted from the top (e.g. if n=10, m=8, hence K n _ m =Ki 0 _ 8 = K 2 ). 

This is the second to the topmost floor or in order words second floor from the top. 

K m = the variable wind pressure coefficient corresponding to the next floor (counted from 

ground floor) to the intermediate floor the building being considered (e.g. if m=8, hence K m = K 8 ). 
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Figure 2.0 (b) : A Section through the Topmost Level of a Multi-Storey Framed 

b.) The forces in columns at the top floor 

The axial forces at the exterior columns (i.e. L= left and R=right) at the top floor are 

Ni = N eLn = -N eRn = - N 4 = r P ln = rP 

The axial forces at the interior columns (i.e. L= left and R=right) at the top floor are 
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N 2 - N intLn - -N intRn - - N 3 - P 2] 


Where, 

Pin - P211 - P 

2 (L^ + O.5L2) 


r = 


...(4) 

•(4a) 

.(4b) 


r = the proporţional ratio between the columns and the center of gravity of the building frame. 
Taking moments about point o;(i. e. ; £ M 0 = 0) 

hh 


Hence, 

P = 


W n --+ pLi - (4 + L 2 )P - (4 + L 2 + LjrP = 0 
h 2 

W n — — P[2rL 1 + (r + 1) L 2 ] = 0 


w n h 2 


•(4c) 


4(2rL^ + (r+l) L 2 ) 

Since Ni = - N 4 = r P 

This is the axial force in the exterior columns at the top most floorof the building. 

Ni = - N 4 = rWnh2 

4(2rL 1 + (r+l) L 2 ) 

and 

N 2 = - N 3 = P 
n 2 = - N 3 = - — W " h2 

4(2rL 1 + (r+l) L 2 ) 

Figure 2.0 (c), shows the horizontal force Hi is determined by taking moments of all internai forces at joint ai 


,(4d) 


Wh/2 


ai 


f 


h/2 


.Hx 


H M al = 0 
h Li 

Hi 2 _Ni y" 0 


N^=rP 


Hx 


0.5Li 


Since Ni = - N 


Figure 2.0 (c) : A Section through the Topmost Fevel of a Multi-Storey Framed Continues 

r W n h 2 


4(2rL 1 + (r+l) L 2 ) 


= Fi 


(4e) 
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Therefore, 


H = 17 rW n h 2 \ Lil 

1 L V4(2rL 1 + (r+l) L 2 )/ hj 

and 

H _ \( r W n L 1 h 

1 l\4(2rL 1 + (r+l) L 2 )/ J 

F 2 = the shear force on tb 
Similarly, to determine H 2 and F 2 ; take sum moment about C 2 . 


\4(2rL 1 + (r+l) L 2 v 

F 2 = the shear force on the beam along the exterior spân of the rigid frame 


W„ h/2 


ai 


h/2 


C"* 

Ni = r P 


F 2 


l 

C 2 


]\f= P 


.H 2 


0.5L/2 J 


.H 2 


Figure 2.0 (d) : A Section through the Topmost Level of a Multi-Storey Framed Continues 
Therefore, 

(i-e.;2M Cz = 0) 

(Hi + H^-Ni (Li+Ş-)-N 2 ^=0 

From this equation, H 2 is the only unknown variable and it can be determined 

H 2 =^(L 1 + if) + N 2 ţf-H, 

Since, 

r W n h 2 . XT _ XT _ W„h 2 


Fi Ni - N 4 4(2rLl+(r+1)L2) > 


N 2 


■N, 


H > = [(; 


rWnLih 




4(2rL 1 + (r+l) L 2 ) 


and 


\4(2rL 1 + (r+l) L 2 )J 

Substituting the values of Ni, N 2 and gives 

H _ r r W n h 2 1 2 / LA r W n h 2 1 L2 _ r rW n L!h 1 
2 [4(2rL 1 + (r+l) L 2 )J h V 1 2/ U(2rL 1 + (r+l) L 2 )J h U(2rL 1 + (r+l) L 2 )J 


H 2 = 

h 2 = 


W n h 


(2rL 1 + (r+l) L 2 ) 
W n h 


4(2rL 1 + (r+l) L 2 ) 


[¥( l . + t) + 7 + t] 

Hi + (1 + r)L 2 ] 


(5) 


(5a) 
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c.) Determination of Axial Force P n _i at the Second to the Last Floor 


Li 

^ 1 kn 

U 

L 






From Figure 3.0: A Section through the Second Floor from the Top of a Multi-Storey Framed Building 




k " K ^" ' (lin + ^y 1 ) + k'w(y + iy 1 ) (yy) + L,P„_, - (L, + L 2 )P„_, - (L, + L 2 + L,)rP„_, = 0 
knWh n _i h n Wh n _ 1 /h n _ 1 \ /h n h n _i\ /h n -i\ T n Ţ 

— 2~ h " + ^—{—) + k w (y — ) + k w (— J + WP-. - O. + «p„-, 

- (L x + L 2 + L 1 )rP n _ 1 = 0 
If h n = h n _! = h ; these equations will give 

k^Wh 2 k^Wh 2 h 2 h 2 

— + k'W— + k'W — + l^P^ - (4 + L 2 )P n _! - (4 + L 2 + LjrP^ = 0 


^(Ş L + k') + LiPn-x - LiPn-i - L 2 P n _ 1 - (L 2 + 2L 1 )rP n _ 1 - 0 
^!(Şl+ k ') - + r) + 2I-LJ = 0 

Wh 2 (^a+k') Wh 2 (3k( 1 +2k') 

n_1 2[L 2 (l+r)+2rL 1 ] 4[L 2 (l+r)+2rL 1 ] 


Finally, 

P„-i = 


Wh 2 (3k{ 1 +2k') = N 
4[L 2 (l+r)+2rL 1 ] 2 


( 6 ) 
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To determine the internai forces at this pre-determined level, take moments along the beam and column at joints ai 
anda 2 respectively. 


d.) At Joint ai 

Along the column 


Since 


Therefore, 


and 

M al 


K 


r k^Lj^Wh 2 
8(2rLi + (r+l) L 2 ). 


;)] 



Hi 


/ rk^Wh V 
y4(2rL 1 + (r + 1) L 2 )J 


M a i 


/ rk^l^Wh \ 
\4C2rL, + (r + 1) L 2 ) / 


h 

2 


Along the beam L,. 

M a i= Fi y 
Fi = N, = -N 4 = 

K: 


r k„ Wh 2 


M 


al 


4(2rL 1 +(r+l) L 2 ) 
rknL x Wh 2 


\8(2rL 1 + (r+l) L 2 ). 


;)] 


Similarly, 

e.) At Joint a 2 

Along the beam L x 


M 


a2 Ll 


- Fi — 
1 2 


M 


a2 Ll 


M 


a2 Ll 


- 1; 


i.) 


rknLiWh 2 1 

8(2rL! + (r+l) L 2 )J 

Along the column h. 


M a2h - H 2 - 


H? = 


knWh 


M 


4(2rL 1 + (r+l) L 2 ) 
k' Wh 


•(7) 


r k„ Wh 2 

L 2 _ 

r 44 Wh 2 

4(21-4 + (r + 1) 4). 

2 

8(2r4 + (r + 1) 4) 


•(7a) 


a2 h 


M 


a2 h 


-H 

-l 


4(2rL! + (r+l) L 2 ) 

knWh 2 


8(2rL 1 + (r+l) L 2 ) 


[rLi + (1 + r)L 2 ] 

l>4 + (1 + r)L 2 ]J 
[rL a + (1 + r)L 2 ]J 


■(7b) 
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ii.) Along the beam L 2 



But 
Ni = - N 4 = 


and 
N 2 = 

N 2 = 


F 2 = Ni + N 2 

rk„Wh 2 


4(2rL 1 + (r+l) L 2 ) 

N 3 = P 


= rP 


■N 3 


k^Wh 2 


4(2rL 1 + (r+l) L 2 ) 


Finally, 


and 


F, = 


rk!,Wh 2 


kI,Wh 2 


4(2rL! + (r + 1) L 2 ) 4(2^ + (r + 1) L 2 ) 


= P(1 + r) 


k„Wh 2 (l + r) 
4(2rL 1 + (r + 1) L 2 ) 



k„Wh 2 (l + r) 
4(2rL 1 + (r + 1) L 2 ) 


h 

2 


' k'Wh z (l+r)L 2 j 
.8(2rLi + (r+l) L 2 )J 


The equilibrium of all moments at joint a 2 is given by 
M a2 h = M a2 L1 + M a2 L2 


M 


a2 Ll 


- 1 


r k^ Li Wh 2 
8(2rL^ + (r+l) L 2 ). 


5 ] 

d aIld M *2L2 


k[ 1 Wh 2 (l+r)L 2 1 

,8(2rL 1 + (r+l) L 2 )J 


M a2 h = 
M a2 h = 


k^Wh 2 


t8(2rL 1 + (r+l) L 2 ) 


[rLi + (1 + r)L 2 ]| = [ 


r k' n L ± Wh 2 


k^Wh 2 (l+r)L 2 


8(2rL^ +(r+l) L 2 ) 8(2rL^ + (r+l) L 2 )J 


k^Wh 2 


l8(2rL 1 + (r+l) L 2 ) 


[rLi + (1 + r)L 2 ]| = 


k^Wh 2 


8(2rL^ + (r+l) L 2 ) 


Hi + (1 + r)L 2 ] 


(7 c) 


Okay 


Generalised Equations for a Determination of Moment and Forces at n th Story Level 


Determination of axial forces at the internai column for the second floor (2 nd ) from ground level (i.e. P 2n ) 


_ [(n-l) + 2£f =1 (n-A)]kWh 2 
2n ~ 4(2rL! + (r+l) L 2 ) 

Total number of storeys n = 10 
Number of storeys j (i.e. n-i) above level i= 
j = n - i = 10-2 = 8 


.( 8 ) 


2 ; 


The last odd number k consider 
k = (ij + 1)= (2x8 +1) = 17 
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[(n- 1) + 2 Zf =1 (n - A)] 

=(n - 1) + 2(n - 3) + 2(n - 5) + 2(n - 7) + 2(n - 9) + 2(n - 11) + 2(n - 13) 

+2(n- 15) + 2(n- 17) 


= 17n-161 


(8a) 


Alternatively, this can be similarly computed using the following procedures 
Sum of all odd numbers = -2j(j+2)-l= -2x8(8+2) -1= -16(10) -1= -161 

Sum of all n = ij+1 = 2x8+1 = 16+1 = 17n (8b) 

Therefore, equation (8) can be re-written as (9) 

n _ [n(ij + 1) - {(2j(j + 2)) + l}]kWh 2 _ [n(2x8 + 1) - {(2x8(8 + 2)) + l)}]kWh 2 
Pj “ 4(2^ + (r + 1) L 2 ) - 4(2rL 1 + (r + 1) L 2 ) 

p _ [17n-16l]kWh 2 _ PiWjh 2 
j 4(2rL 1 + (r+l) L 2 ) 4(2rL 1 + (r+l) L 2 ) 

If P n — 1 + 2(nj) 

and since j = n — 1; 

Therefore, let 

p n = 1 + 2(nj) or p n = 1 + 2n(n - 1) or p n = 2n 2 - 2n + 1 

/n + 3\ /2n - 3\ 

w n = (W min + (—) K max )b = 0.1042bW bas ; V,n > 3 

- h 2 

a " ~~ 4[2rL!+(l+r)L 2 ] 


(9a) 

(9b) 


From Figure 4.0a, the axial load, Pi„at the left and right exterior columns is given by 


Pin 


r[2n 2 -2n+l][(W m i n +(^)K max )b]h 2 
4[2rL 1 + (r+l) L 2 ] 


(9c) 


and the axial load,P 2 „at the left and right interior columns is given by 


(1 + 2n(n - 1)} [(W min + (^) K max )b] h 


P,n = 


P,n = 


_ [2n 2 -2n+l][(W min + (H±l) Kmax )b]h 2 


4(2rLi + (r + 1) L 2 ) 


4[2rLi + (r+l) L 2 ] 


(9d) 


From Figure 4.0 (a), all these values of Pi n and P 2n (in 9c & 9d) can be re-written for the n th storey counted from the 
topmost floor as equations (9e and 9f) 


p _ r p n w n h 2 

ln 4(2rLi+(r+l) L 2 ) 
p _ p n W n h 2 
2n 4(2rL!+(r+l)L 2 ) 


f a nPnW n 

a n p n w n 


(9e) 

(9f) 
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Shearing forces H ln and H 2n along the columns are 

(2n-l)rW n h 2 _ r (2n-l) anWnLi 


H 1n = 


H? n = 


4(2rL 1 +(r+l) L 2 ) 
(2n-l)W n h 2 [rL 1 + (l+r)L 2 ] 
4(2rL 1 + (r+l) L 2 ) 


a n W n (2n-l)[rL 1 +(l+r)L 2 ] 


Shearing forcesF ln and F 2n along the beams are 


Fm = 


F 2n = 


(n-l) rW n h 2 
4(2rL 1 + (r+l) L 2 ) 
(n-l)(r+l) W n h 2 
4(2rL 1 + (r+l) L 2 ) 


(n - l)4ra n W n 
(n — l)(r + l)4a n W n 


(10a) 

(10b) 


(lla) 

(llb) 


Bending moments along the beams and columns shown in Figure 4.0(b), 

The moment at underside (i.e. bottom) of an exterior rigid joint of a column for a particular storey level, n is M Pnc 

Mfnt = Hin î =n sn^^ (,2a) 

The moment at topside (i.e. top) of an exterior rigid joint of a column for a particular storey level, n is M Pnc 

Mine = H,„-4 = 'C—y-W-L, (,2b) 

The moment along the exterior rigid joint of a beam for a particular storey level, n isM lnb 

M lnb = F ln ± = 2r(n - l)a n W n L 1 (12c) 

But, when n=l,the moment along the exterior beam at the topmost floor M llb is given by 

M llb = F 11 y = \ anWnL-L (12c) 

The moment at underside (i.e. bottom) of an interior rigid joint of a column for a particular storey level, n is M 2nc 
M 2 B „ c = H 2n \ = ^W n (2n-l)[rL 1+ (l + r)L 2 ] (12d) 

The moment at topside (i.e. top) of an exterior rigid ioint of a column for a particular storey level, n is M J nr 

m T h h a n W n (2n-l)[rL 1 +(l+r)L 2 ] 

M 2nc = H 2n -i - = (12e) 

The moment along the interior rigid joint of a beam for a particular storey level, n isM 2nb 

M 2 nb = M 2 R nb = F 2n = 2(n - l)(r + l)a n W n L 2 (12d) 


Analytical Examples 

A 14m three bays (F!=6m, L 2 =4m) framed medium rise 48m (h=3m), 16 storey building is 3 Om long with 6m 
columns or frames spacing (i.e.six transverse braced frames). Basic wind speed for Maiduguri where the building is 
sited is W bas =47 ms" 1 ; altitude factor S a = 1.34 for an elevation of 340m above the sea level and the following 
relevant factors: permanent seasonal S p = 1.0; directionality S d = 0.99, dynamic K b = 1.0;dynamic augmentation C r 
=0.04; size effect C a = 0.83; externai surface pressure C pe = 1.1; and terrain factor varying between 1.362 at lOm 
and 1.900 at 48m above ground level. The conversion factor of a triangular pressure to evenly distribution pressure 
is 0.625. Using the appropriate equations and formulae determine the axial, shearing forces and their corresponding 
bending moments on the columns and beams; at the sixth floor of a rigid unbraced frame subjected mainly to 
aerodynamic loadings. 


Solution 

From equation 2f, the equivalent uniformly distributed load from a W bas = 47 m/s applied at each floor level is 


w„ = W m „ = (W m „, + g K lli;i , ,l> = 0.1042bW„„ gg) ; V, n > 3 
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W n = 0.1042x6x47x 


2x16 - 3 
.16-3 


) 


= 65.53 kN 


Fin F2 " F ln 



Figure 4.0 (a) : The Axial and Shearing forces along the Columns and 



Figure 4.0 (b) : The Bending Moments on the Columns and Beams 


Therefore, column forces and moment are determined as follows 

a.) The Exterior, P ln and Interior, P 2n Column Axial Forces 

From equations 9(a and b), the non-dimensional storey number (3 n and variable dimensions, a n influence 
coefficients are determined. 

p n = 2n 2 — 2n + 1 = 2xl6 2 x2 - 16 + 1 = 481 
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h 2 


a n = 


4[2rL x + (1 + r)L 2 ] 4[2x4x6 + (1 + 4)4] 


= 0.03309 


From equations 9(c, d and f) 


r P n W n h 2 

ln 4(2rL x + (r + 1)L 2 ) 

PnW n h 2 

4(2rL! + (r + 1) L 2 ) 


= r a n (B n W n = 4x0.03309x481x65.53 = 4172 kN 
a n p n W n = 0.03309x481x65.53 = 1043 kN 


b.) 


H 1n = 


The Exterior, Hi„ and Interior, H 2n Shearing Forces along the Columns 

From equations 10(a and b), 


(2n-l) r W n h 2 
4(2rL 1 + (r+l) L 2 ) 


r (2n-l) anWnLi _ 4x0. 03309x65. 53x6x(2xl6-l) 


= 538 kN 


H? n = 


(2n - l^h 2 ^ + (1 + r)L 2 ] = a n W n (2n - ljţrL, + (1 + r)L 2 ] 
4(2rL x + (r + 1) L 2 ) h 


TT 4x0, 03309x65. 53(2xl6-l)x[4x6+(l+4)4] 8.674x31x44 

H 2n = — 1 — — = = 986 kN 


c.) The Exterior, F in and Interior ,F 2n Shearing Forces along the Beams 

From equations 1 l(a and b), 


Shearing forcesF ln and F 2n along the beams are 

,W n = (16 
(n — l)(r + l)W n h 2 
4(2rLi + (r + 1)L 2 ) 


Fln = 4 ( 2 rLi+(r+i) 1 L 2 ) = ( n “ !) 4ra n w n = (16 - I)x4x4x0.03309x65.53 = 520.41 kN 


F?n = 


= (n — l)(r + l)4a n W n 


F 2n = (16 - 1) (4 + I)x4x0.03309x65.53 = 650.52 kN 


d.) Bending Moments along the Beams and Columns 

From equations 12(a - d) 

The moment at underside (i.e. 16 th Floor) of an exterior rigid joint of a column for a particular storey le vel, n is M 

Mfnc = H ln ^ = r a " W " L ^ 2n -^ = 4x0.03309x65.53x6x(2xl6— 1) = QQJ RNm; ( . Hi = 533 ^ 5) 


The moment at topside (i.e.!5 th Floor) of an exterior rigid joint of a column for a particular storey level, n is Mf nc 


Mine 


= H ln -i- = 


r a n W n L 1 (2n-l) _ 4x0. 03309x65. 53x6x(2xl5-l) 


= 754kNm 


The moment along the exterior rigid joint of a beamfor a particular storey level, n isM lnb 

M lnb = F ln F = 2ra n W„L 1 (n - 1) 

M lnb = 2x4x0. 03309x65. 53x6x(16 - 1) = 1561 kNm (i. e. 520.41 x3 ) 


But, when n=l,the moment along the exterio r beam at the topmost floor M llb is given by 
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Mnb = Fu Y = T - anWnLi = Jx0.03309x65.53x6 = 26 kNm 

The moment at underside (i.e. 16 th Floor) of an interior rigid joint of a column for a particular storey level, n is Mf nc 

M 2 B nc = H 2n t = a n W n (2n-l)[rL 1 + (l + r)L 2 ] = 0.03309x65.53x(31)x[24 + 20] = ^ ^ ( . 9 g 6xl5) 

The moment at topside (i.e. 15 th Floor) of an exterior rigid ioint of a column for a particular storey level, n is Mj nr 

wT _ TÎ h_ a n W n (2n-l)[rL 1 + (l + r)L 2 ] 

M 2nc — H 2n-1 2 

AyrT 0.03309x65. 53(29)[44] 100 oi m r* om i ca 

M 2nc = = 1383kNm ( 1 . e. 922x1.5) 

The moment along the interior rigid joint of a beam for a particular storey level, n isM 2nb 

M 2 nb = M 2 R nb = F 2n F - 2(n - 1) (r + l)a n W n L 2 
M 2nb = 2x(15)(5)x0.03309x65.53x4 = 1301 kNm (i. e. 650.52x2) 

e.) Equilibrium of Joints 

From Figure 4b; in complying with the requirements of the principie of minimum potenţial energy at all joints, the 
joint moment equilibrium is given as folio ws. 

For the exterior columns, 

Mine + Mfnc + M lnb = 807 + 754 - 1561 = 0 
Similarly, for the interior columns 

M lnb + M 2nb - Mţ bc - M 2nc = 1561 + 1301 - 1383 - 1479 = 2862 - 2862 = 0 
Similar equilibrium is obtained from sum horizontal and vertical force for the substitute frames. 

DISCUSSION OF RESULTS 

The slab is considered to be an infinitely rigid structural element along the plane of the frame; its rigidity assists in 
the distribution of the horizontal load from wind or seismic loadings to be concentrated not within the inter storey 
column height but at the levels of the slab which acts as a non-deformable tie to unite all the frames together with 
any available primary and secondary beams. Although, for most practicai situations, the need for beams and 
columns are usually dictated by the architectural and technological considerations, but for structural purposes, the 
sizes of reinforced concrete beams in millimeters are to be determined by equation 13. 



where; M 0 is the bending moment at the beam joint for wind or seismic loadings or mid spân bending moment for a 
simply supported beam from vertically acting dead and imposed loadings; b > ^ is the concrete web, F is the beam 

effective spân and f ck characteristic strength of the concrete. Depending on the percentage of reinforcement required, 
the required area of concrete column may be determined by equations 14(a and b) 

A c = —for N < lOOOkN ... 14a 

fck 

A c = ^31 for N > lOOOkN but Ac > 122,500 nun 2 ... 14b 

fck 


From the bending moment of 1479 kNm obtained at the 16 th floor counted from the top, assuming a beam web b w 
=300mm and grade 25 concrete, from equation (13), a 650mm overall depth h of the beam may be required. 
Similarly, since the exterior column gave an axial force of 4172 kN, eq. (14b) shows the area of appropriate concrete 
concrete is 375480mm 2 . This could squared column of 625mm or 500x750mm respectively. This approach and 
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values offer useful preliminary criteria for preliminary sizing of beams and columns for rigid frames subjected to 
aerodynamic loadings. 


REFERENCES 

ASCE Standard, ASCE7-95(1995). Minimum Design Loads for Buildings and Other Structures, American 
Society of Civil Engineers, New York. 

British Standard, Loadings for Buildings (B.S. 6399) - Part 2: Code of Practice for Wind Loads (as Amended) 
(2004), Building and Civil Engineering Sector Board, United Kingdom. 

Grider, A., Ramirez, J.A. and Yun, Y.M. (1999). “Structural Engineering Handbook”. Structural Concrete 
Design.pdf. Fourth Edition. Edited by Chen Wai-Fah. Boca Roton: CRC Press LLC. 

Joseph, A.M., William PF (2006). Database-Assisted Design for Wind: Concepts, Software, and Examples for Rigid 
and Flexible Buildings, NIST Building Science Series 180. Special Case Wind Design, Main Fritz _DAD_BSS180 
Pdf. 

Kijewski, T, Kareem A (2001). Dynamic Wind Effects: A Comparative Study of Provision in Codes and Standards 
with Wind Tunnel Data, Tracy.L.Kijewski.l@nd.edu , Notre Dame. 

Mendis, P, Ngo T, Haritos N, Hira A, Samali B, Cheung J (2007). Wind Loading on Tall Buildings, EJSE 

Mosley, W.H., Bungey, I.H.and Hulse, B. (2007). “Reinforced Concrete Design”. Sixth Edition, Published by 
Palgrave. New York. 

Onundi, L.O., Oumarou, M.B. and Philip Otuagbu. (2009). Classification of Nigeria Wind Isopleths for Structural 
Purposes. Continental J. Engineering Sciences. Wilolud Journals. 4:48 - 55. Available online at 
http : // www . wiloludi ournal . corn . 

Onundi ,L. O., Elinwa A. U. and Matawal D. S.(2012). Development of a wind tunnel model compared with results 
of method of iniţial parameters and BS6399 models for dynamic analysis of multi-storey building subjected to 
aerodynamic loadings. Journal of Civil Engineering and Construction Technology 3(5), 159-178. Available online 
at http://www.academiciournals.org JCECT DOI: 1 0.5 897/JCECT 11.064 .ISSN 2141-2634 ©2012 Academic 
Journals 


Onundi, L.O. (2010). The Impact of Climate Change on Sustainable Infrastructural Development - A Case Study of 
the Appropriate Wind Speed and other Measures Required for Design of Tall Structures In Nigeria. A Paper 
Presented to the Council for Regulation of Engineering in Nigeria 2010 Assembly (2010 COREN Assembly). 
Abuja. 

Oyenuga, V. O. (2011). “Simplified Resinforced Concrete Design (A Consultant / Computer Based Approach). 
Second Edition. Published by Asros Limited. ISBN : 978-36217-7-7.Lagos. 


Smith, B.S. and Coull, A. (1991). “Tall Building Structures: Analysis and Design”, Published by Wiley Interscience 
Publications. John Wiley and Son Inc. ISBN : 0-471-51237-0. New York. 

Special Issue: Loading on Structures. pp. 1-14. 

Soboyejo, A.B.O. (1971). Distribution of Extreme Winds in Nigeria, Niger. Eng., 7(3): 21-34. 


Taranath, SB (2005). Wind and Earthquake Resistant, Structural Design and Analysis, Marcel Dekker, New York. 

|@ 0 ~ 1 

AII rights reserved 

This work by Wilolud Journals is licensed under a Creative Commons Attribution 3.0 Unported License 

26 



